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Abstract: The non-rectangular hyperbola (NRH) equation is the most popular method that plots the 
photosynthetic light-response (PLR) curve and helps to identify plant photosynthetic capability. However, 
the PLR curve can't be plotted well by the NRH equation at different plant growth phases due to the 
variations of plant development. Recently, plant physiological parameters have been considered into the 
NRH equation to establish the modified NRH equation, but plant height (H), an important parameter in 
plant growth phases, is not taken into account. In this study, H was incorporated into the NRH equation 
to establish the modified NRH equation, which could be used to estimate photosynthetic capability of 
herbage at different growth phases. To explore photosynthetic capability of herbage, we selected the 
dominant herbage species Potentilla anserina L. and Elymus nutans Griseb. in the Heihe River Basin, 
Northwest China as the research materials. Totally, twenty-four PLR curves and H at different growth 
phases were measured during the growing season in 2016. Results showed that the maximum net 
photosynthetic rate and the initial slope of PLR curve linearly increased with H. The modified NRH 
equation, which is established by introducing H and an H-based adjustment factor into the NRH equation, 
described better the PLR curves of P. anserina and E. nutans than the original ones. The results may 
provide an effective method to estimate the net primary productivity of grasslands in the study area. 
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1 Introduction 


The photosynthetic light-response (PLR) curve describes the relationship between the available 
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photosynthetic photon flux density (PPFD) and the leaf net photosynthetic rate (PN) (Amanullah, 
2015; Noda et al., 2015). The PLR curve is an analytical tool that helps identify the important 
feature of plant, such as photosynthetic capability. Up to date, numerous models or equation have 
been constructed to describe it. For example, quadratic polynomial regression model, subsection 
model, rectangular hyperbola model, non-rectangular hyperbola (NRH) equation, exponential 
model, modified rectangular hyperbola model and modified exponential model are used to 
simulate the PLR curve (Thornley, 1976; Farquhar et al., 1980; Marshall and Biscoe, 1980; Kull 
and Kruijt, 1998; Zufferey et al., 2000; Ye, 2007; Ye et al., 2013; XU et al., 2014; Huntingford et 
al., 2015). Among these methods, the NRH equation is the most popular one (Prieto et al., 2010; 
Calama et al., 2013; XU et al., 2014; Mayoral et al., 2015; Zvalinsky and Tishchenko, 2016). 

The photosynthetic capability is influenced by PPFD and physiological and morphological 
features of leaf (Åström et al., 2015; Catoni et al., 2015; Sakoda et al., 2016; Way et al., 2017). 
The physiological and morphological features of leaf generally include chlorophyll (Chl) content 
(Walker et al., 1983; Keutgen et al., 1997; Zhu et al., 2012; Gitelson et al., 2015; Khoshbakht et 
al., 2015; Novichonok et al., 2015; Sawicki et al., 2017; Zhu et al., 2017), water status (Quero et al., 
2008; Akhkha, 2010), nitrogen (N) content (Catoni et al., 2015; Onoda et al., 2017), age 
(Chiarawipa et al., 2012), specific leaf mass (Puglielli et al., 2015, 2016; Xiong et al., 2016) and 
CO: concentration (Urban, 2003; Wang et al., 2015, 2017; Smith et al., 2016; Hao et al., 2017). 
Therefore, modified PLR curves are various due to different environmental conditions of each 
plant species (Iacono and Sommer, 2010; Lachapelle and Shipley, 2012; Legner et al., 2013; Mu 
and Fang, 2015; Malerba et al., 2017). 

Recently, some researchers are beginning to consider leaf factor into the PLR curve. Ohsumi et 
al. (2007) incorporated leaf N content into the equation to calculate the photosynthetic rate of rice 
leaves. Marino et al. (2010) and Lachapelle and Shipley (2012) established an allometric equation 
to estimate Py with specific leaf mass and N content. Calama et al. (2013) incorporated leaf 
temperature, soil moisture, and leaf type into the NRH equation and predicted the Py of pine 
needles. Xu et al. (2014) established a modified NRH equation by incorporating Chl and leaf N 
contents to estimate the soil plant analysis development (SPAD) value. However, the 
determinations of these physiological factors are always destructive and time consuming. We 
attempted to find a morphological factor to modify the NRH equation to save time and resources. 
Plant height (H) is a popular indicator to reflect the growth of plants (Ansquer et al., 2009; 
Muharam et al., 2014; Teh et al., 2015; Xing et al., 2015; Cao et al., 2017; Gleason et al., 2017; 
Querry and Harper, 2017). However, H is rarely considered in the NRH equation to estimate Py 
for grass species. 

Potentilla anserina L. and Elymus nutans Griseb. are dominant species of herbages in the 
Heihe River Basin, Northwest China. The values of net primary productivity (NPP) of the species 
have received great attention due to development of the livestock industry. Moreover, NPP is 
closely related to photosynthetic capacity (Malhi et al., 2017). The photosynthetic characteristics 
of grass are referred to as specific leaf mass, water stress, light condition, leaf N content, leaf 
position and H. Among these characteristics, H is easily and non-destructively measured and it 
can represent leaf photosynthetic capacity during the growing season. Thus, our objective is to 
incorporate H into the NRH equation to establish the modified NRH equation, which can be used 
to describe the PLR curves of herbage at different growth phases. The results may provide an 
effective method to estimate the NPP of grasslands in the Heihe River Basin, Northwest China. 


2 Materials and methods 


2.1 Study area 


The Heihe River, the second largest inland river in the arid region of Northwest China, originated 
from the Qilian Mountains. Our study area is located in the upper reaches of the Heihe River 
Basin (38°23'56"—38°26'47"N, 99°53’57"—99°57'10"E), covering an area of approximately 12.8 
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km? with an elevation range of 2660-4419 m a.s.1. Its climate is characterized by short and warm 
summer and long and cold winter. The annual mean temperature is from —0.6°C to 2.0°C, and the 
mean annual precipitation is 435.5 mm. Temperature and precipitation vary with the elevation, i.e., 
temperature decreases and precipitation increases with increasing elevation. The variation causes 
a distinct vertical distribution of plant communities from low to high elevation as follows: desert 
steppe, forest steppe, subalpine shrub and subalpine meadow (Gao et al., 2016). 


2.2 Experiment and measurement 


The study area belongs to subalpine meadow. The dominant species are P. anserina and E. nutans. 
Two 8 mx8 m plots, one for P. anserina and the other for E. nutans, were selected before the 
growing season in 2016. Each plot was divided into 4 subplots (4 mx4 m). A 10 m-wide buffer 
was designed between two plots as protection from the disturbance of edge effect. The PLR 
curves of leaf were measured using the Li-6400 portable photosynthesis system (LI-COR Inc., 
Lincoln, NE, USA) at the early growth phase (15 June), the middle growth phase (3 and 22 July) 
and the late growth phase (15 August) in 2016. In each measurement, the healthy and fully 
expanded leaves at the middle-upper part of the plant were selected. The selected leaves were 
measured after 20 min of illumination at a PPFD of 2000 umol/(m*s). Px was obtained at 14 
levels of PPFD (2000, 1800, 1500, 1200, 1000, 800, 600, 400, 200, 150, 100, 50, 20 and 0 
umol/(m?+s)) under natural conditions. One subplot was selected at a time to measure the PLR 
curve for each species. Three PLR curves can be drawn according to the corresponding PPFD of 
Py. After measuring the PLR curves, H was measured, and 9 replicates were obtained in each 
subplot. The PLR curve and H were measured 4 times during the growing season, resulting in 24 
PLR curves. 


2.3 Development of the modified NRH equation 
Unmodified NRH equation is expressed as follows (Marshall and Biscoe, 1980; Thornley, 1998): 


al+P. (al+P > _ A@a IP. 


max E max ) max 
P (I) = N Ti m L Rp (1) 


where a is the initial slope of the PLR curve (°); I is the PPFD (umol/(m7ss)); 6 is the convexity of 
the PLR curve; PNmax is the maximum net photosynthetic rate (umol/(m?+s)); and Rp is the dark 
respiration (umol/(m7+s)). 

Parameters of a, 0, PNmax and Rp were calculated by non-linear least-square fitting. And we 
developed the modified NRH equation by incorporating H into the equation: 


al + Pina — COL + Poma) — AOGTP 
Py (1) = BH Nmax ( aa! Nmax Ry. (2) 


where B is a coefficient; and H is the height of species (cm). The mean values of a and PNmax in 
the modified NRH equation were used to define B. 


2.4 Statistical analysis 


We used Sigmaplot (Version 12.5, Systat Software Inc., USA) statistical packages to analyze 
measurements. The parameters in Equations 1 and 2 were determined by non-linear least-square 
fitting, which were realized by software 1“Opt (First Optimization; 7D-Soft High Technology Inc., 
China) using the Levenberg-Marquardt (LM) algorithm. The average absolute error (AE) and the 
root mean square error (RMSE) were used to assess the performance of the unmodified NRH and 
the modified NRH equations, and they were calculated as follows: 


1 n 
AE=— D , 6 


leon 
RMSE = E > Pi = Fhovs,i y > (4) 


where PNeaii is the calculated Pyn by Equations 1 and 2 (umol/(m*+s)); PNobs,i is the observed Py 
(umol/(m*-s)); and n is the total number of data (n=336). 


Peat = Pyobs.i 
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The performance of the modified NRH was based on 4 prioritization schemes, that is, Scheme 
1: Inputting Py and PPFD into the modified NRH equation; Scheme 2: Inputting the a limited 
from 0.03 to 0.08 into the modified NRH equation; Scheme 3: Inputting the maximum values of 
Pymax and a into the modified NRH equation; Scheme 4: Inputting the mean values of Pymax and a 
into the modified NRH equation. According to the measured and estimated photosynthetic 
parameters, we selected the best scheme. 


3 Results 
3.1 PLR curves and H at different growth phases 


The PLR curves of P. anserina and E. nutans were shown in Figure 1. The mean Py values were 
significantly higher in July and August than those in June. The difference in Py at various growth 
phases was small when the PPFD was below 500 umol photon/(m?-s), and the difference in Py 
increased when the PPFD was above 500 umol photon/(m?+s). Figure 1 showed that the maximum 
Py values of P. anserina on 15 June and 15 August were 11.54 and 26.71 umol CO2/(m*ss), and 
the maximum Py values of E. nutans on 15 June and 15 August were 12.17 and 18.03 umol 
CO,/(m*+s), respectively. The variation in H with time of the two species was shown in Figure 2. 
The maximum and minimum H values for P. anserina were 16.62 and 6.13 cm and those for E. 
nutans were 68.99 and 36.09 cm, respectively. Moreover, as shown below, high H corresponded 
with high Py rate. 
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Fig. 1 Net photosynthetic rate (Pn) response to available photosynthetically photon flux density (PPFD) for the 
leaves at various growth phases in 2016. Bars mean standard errors; n=3. 
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Fig. 2 Variation of plant height (H) values of the two species at various growth phases in 2016. Bars mean 


standard errors; n=9. 
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3.2 Relationship between parameters of the NRH equation and H 


For each PLR curve, the NRH equation fitted very well, because photoinhibition of plant species 
was not found (Fig. 1). The parameters of the NRH equation varied with the changes of H at 
different growth phases. By analyzing the relationship between photosynthetic parameters in 
Equation 1 and H, we found that two photosynthetic parameters (i.e., æ and PNmax) increased 
linearly with H. And the linear relationships were significant at the 0.99 confidence level. 
Furthermore, the parameters of 0 and Rp were not significantly related to H (Fig. 3). 
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Fig. 3 Relationships between plant height (H) and photosynthetic parameters (a, PNmax, 0 and Rp) in the NRH 
(non-rectangular hyperbola) equation. “ indicates the significance at P<0.01 level. a, initial slope of the PLR 
curve; PNmax, maximum net photosynthetic rate; 0, convexity of the PLR curve; Rp, dark respiration. 


3.3 Performance of the unmodified and modified NRH equations 


Comparing the values of observed and estimated photosynthetic parameters, we found that 
Scheme 4 (i.e., inputting the mean values of PNmax and a into the modified NRH) was the best 
because the values of parameters really reflected the growing process of the two species (Table 1). 
Moreover, the mean values of PNmax and a are similar to those of the modified NRH equation in 
Scheme 4 (Table 2). Therefore, mean values of PNmax and a were selected to run the modified 
PLR curves. For P. anserina, the values of 0.0488° (a) and 23.2359 umol/(m*+s) (PNmax) Were set as 
the mean values for the 12 PLR curves, while 0.0572° (a) and 19.1421 umol/(m?+s) (PNmax) were set 
for E. nutans (Table 2). Values of 0, Rp, and B in Scheme 4 for P. anserine were determined as 
0.7089, 0.8436 umol/(m*-s) and 0.0767, respectively, while for E. nutans they were 0.4749, 1.2876 
umol/(m*+s) and 0.0172, respectively (Table 1). Moreover, we found that the performance of the 
modified NRH equation performed better than that of the unmodified NRH equation with a higher 
R? value (Fig. 4). 
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Table 1 Values of parameters from different prioritization schemes by the 1“Opt (First Optimization) software 
with Levenberg-Marquardt algorithm for P. anserina and E. nutans 


Modified NRH equation 


Species Parameter Unmodified 
NRH Scheme 1 Scheme 2 Scheme 3 Scheme 4 
0 0.7453 0.8643 0.8643 0.8150 0.7089 
Rp (umol/(m?.s)) 1.7223 0.2926 0.2926 0.5287 0.8436 
Pymax (umol/(m?-s)) 23.1512 2.7814 18.6565 31.6580 23.2360 
P. a (°) 0.0483 0.0046 0.0311 0.0580 0.0488 
anserina RMSE 3.5513 1.3235 1.3235 1.3411 1.4148 
AE 2.3748 1.1063 1.1063 1.1178 1.1733 
R 0.8377 0.9775 0.9775 0.9769 0.9742 
p - 0.6100 0.0909 0.0546 0.0767 
0 0.5123 0.6488 0.6488 0.1044 0.4749 
Rp (umol/(m°.s)) 1.8234 0.9835 0.9835 1.6033 1.2876 
Pymax(umol/(m*-s)) 19.0431 0.3028 11.8002 21.3590 19.1421 
a (°) 0.0562 0.0008 0.0300 0.0782 0.0572 
eee RMSE 1.7629 1.0400 1.0400 1.1179 1.0587 
AE 1.3273 0.8100 0.8100 0.8281 0.8058 
R 0.9269 0.9746 0.9746 0.9706 0.9737 
p - 1.0501 0.0269 0.0162 0.0172 
Note: 0, convexity of the PLR curve; Rp, dark respiration; Pnmax, maximum net photosynthetic rate; a, initial slope of the PLR curve; 
RMSE, root mean square error; AE, absolute error; R?, coefficient of determination; B, coefficient. "-" means no value. NRH, 


non-rectangular hyperbola. The four schemes are presented in Section 2.4. 


Table 2 Parameters variations in the means, unmodified and modified equations for P. anserina and E. nutans 


Species Parameter Mean Unmodified NRH equation Modified NRH equation 
a (°) 0.0488+0.0071 0.0483 0.0488 
Pymax (umol/(m*-s)) 23.2359+6.1576 23.1512 23.2359 
P. anserina 
0 0.7384+0.0914 0.7452 0.7089 
Rp (umol/(m?.s)) 1.6535+0.5493 1.7223 0.8436 
a (°) 0.0574+0.0129 0.0562 0.0574 
Pyimax (umol/(m?-s)) 9.1421+2.2305 19.0431 19.1421 
E. nutans 
0 0.4974+0.1220 0.5123 0.4750 
Rp (umol/(m°.s)) 1.8458+0.4643 1.8234 1.2876 


Note: a, initial slope of the PLR curve; PNmax, maximum net photosynthetic rate; 0, convexity of the PLR curve; Rp, dark respiration; 
Mean+SD; n=12. 


4 Discussion 


The NRH equation describes the relationship between Py and PPFD (Mayoral et al., 2015; 
Zvalinsky and Tishchenko, 2016), which is popular for expressing the PLR curves. However, it is 
laborious and time consuming to estimate photosynthetic capability (Xu et al., 2014). A lot of 
studies have depicted photosynthetic characteristics of species on the diurnal scale (Belshe et al., 
2008; Orujova et al., 2018), few studies reported photosynthetic characteristics on the seasonal 
scale (Proietti and Famiani, 2002; Gao et al., 2015). However, estimate of Pn depended on limited 
measurements cannot represent the photosynthetic capability during plant growth phases. 
Therefore, studies have tried to estimate the trend of photosynthetic capability by measuring 
biological factors and combining with parameters in the NRH equation (Ohsumi et al., 2007; 
Lachapelle and Shipley, 2012). In this way, Pn could be interpolated during growth season, and 
the accuracy of Px estimated would be improved. 
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Fig. 4 Relationships between observed (PN-observed) and calculated Py (PN-caiculatea) values. NRH, non-rectangular 
hyperbola. 


Four parameters (i.e., æ, 0, PNmax and Rp) in the NRH equation are generally obtained by fitting 
of measured PLR curve. Several group values of a, 0, PNmax and Rp were derived according to 
measuring methods. If the correlation between parameters and biological factors was low, the 
range of parameters could be selected, otherwise, certain values were selected (Marino et al., 
2010; Xu et al., 2014). For example, Marino et al. (2010) selected the certain values of PNmax and 
Rp to input the NRH equation, but Xu et al. (2014) selected the ranges of 0 and Rp to input the 
NRH equation. We selected the ranges of @ and Rp, and certain values of a and Pnmax in the study. 
Regarding a and Pymax, Xu et al. (2014) directly selected the maximum values of æ and Pnmax to 
run the SPAD-NRH equation. After analyzing the 4 schemes, we chose the mean values of a and 
Pymax to the modified PLR curves, which really reflected the growing process of the two species. 

Easily measured factors such as H, leaf Chl content, N content, leaf water content and leaf 
special mass are highly related to photosynthetic capability (Akhkha, 2010; Scartazza et al., 2016; 
Bielczynski et al., 2017; Husen, 2017). Therefore, leaf Chl content, N content, leaf water content 
and specific leaf mass have been introduced into the NRH equation to improve the accuracy of 
estimation (Marino et al., 2010; Xu et al., 2014). Comparatively, H can be measured more easily 
and non-destructively. In our study, we attempt to introduce H into the NRH equation to develop a 
new equation, that is, the modified NRH equation. The modified NRH equation has two 
advantages. On one hand, it had a better performance compared to the unmodified NRH equation. 
On the other hand, it overcame the shortcoming of the unmodified PLR curves, which is time and 
labor consuming. 

The performance of the modified NRH equation is also better than that of the SPAD-NRH 
equation (R?=0.9686; n=10), in which the parameters are calculated using the Chl and N contents 
of leaf (Xu et al., 2014). Thus, H can estimate the Py of plants better because Py can be simulated 
throughout growing season. Furthermore, this method is easy to apply because its parameters can 
be easily obtained. Consequently, it is more useful for describing photosynthetic responses under 
different growing conditions. 


5 Conclusions 


The PLR curves and parameters of herbage in the NRH equation changed with different values of 
H. In the NRH equation, a and Pnmax increased linearly with H. The modified NRH equation is 
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established by multiplying the linear adjustment factor and H (PXH) in the NRH equation. The 
modified NRH equation has an advantage over the unmodified NRH equation in terms of 
accuracy. The performance of the modified NRH equation should be further validated by 
measuring more PLR curves using different species and growing phases. The study of combining 
the corresponding environmental factors into the NRH equation would extrapolate to other 
species or locations. 
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